Ferroelectricity induced by cooperative orbital ordering and Peierls instability. 
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A general mechanism by which orbital ordering, coupled to Peierls-hke lattice distortions, can 
induce an electronic switchable polarization is discussed within a model Hamiltonian approach in 
the context of the modern theory of polarization. By using a Berry-phase approach, a clear picture 
emerges in terms of Wannier-function centers and orbital occupancies. The proposed mechanism 
may apply to oxide spinels whose electronic structure has effective one- dimensional character, such 
as CdV2 04, recently proposed to display multiferroic behaviour. 
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Introduction. Since the revival of interest in the mag- 
netoelectric effect and in the multiferroic state of con- 
densed matter the possible electronic, as opposed 
to purely ionic, origin of ferroelectric polarization has at- 
tracted, and still attracts, increasing attention^, A 
significant contribution to the field came from the rig- 
orous definition of macroscopic bulk polarization P in 
crystal lattices, as given in the early 90's, which eventu- 
ally led to the formulation of what is nowadays known 
as the modern theory of polarization0-Q. According 
to this theory, polarization changes can be formally ex- 
pressed in terms of the Berry phase (BP) of the electronic 
wavefunctions, which correctly captures the contribution 
to P coming from mobile electrons in crystals with co- 
valent character - what is generally referred to as the 
electronic contribution to polarization. Several mecha- 
nisms responsible for electronic ferroelectricity has been 
devised since, involving spin[i charge [llj, [l^ and 



orbital|13L Il4| degrees of freedom of electrons. 



Due to its potential straight application in magneto- 
electric devices, the first class of mechanisms, that ap- 
plies directly to magnetically induced ferroelectrics, has 
been object of intense research^, 0, [l^. Nonetheless, 
our understanding of general mechanisms leading to elec- 
tronically induced ferroelectricity cannot neglect the im- 
plications of charge- (CO) or orbital-ordering (00) phe- 
nomena. The essential mechanism by which a charge 
ordering can lead to a bulk polarization is fairly well 
understood 0, |T3|- Ferroelectricity, in fact, may arise 
due to a simultaneous inequivalency of both sites and 
bonds after charge ordering, a situation often found in 
transition metal compounds containing transition metal 
ions with different valence. Bond inequivalency may 
be structural, as for example in quasi one-dimensional 
organic materials (such as (TMTTF)2X[il| or TTF- 
CA[i6|). or may have a magnetic origin, as in the case 
of CaaCoMnOe^ or i?Ni03[ii|, where ferromagnetic 
and antiferromagnetic bonds differentiate through sym- 
metric exchange striction. The simultaneous presence of 
site-centered and bond-centered CO has been also pro- 
posed to induce ferroelectricity in half-doped manganites 



such as Pro.sCao.sMnOa or Lao.sCao.sMnOaflll. In this 
case, however, we have suggested that orbital degrees 
of freedom may play a relevant role in establishing the 
ferroelectric state [19|. The interplay of electron-electron 
(e — e) and electron-lattice (e — I) interaction, in partic- 
ular the coupling with the buckling mode associated to 
rotation of the oxygen cages ubiquitous in perovskyte ox- 
ides, causes an effective dimerization via an orbital stric- 
tion mechanism, also responsible for the predicted 00 
pattern Quite surprisingly, 00 combined with 

the forementioned effective dimerization appears to be 
the main source of an electronically induced polarization, 
since CO is strongly inihibited in the dimerized statefljj. 

In our opinion, these findings call for further research 
and suggest that the role of 00 in establishing or con- 
tributing to ferroelectric polarization deserves a deeper 
understanding. In order to provide a clear picture of the 
possible general mechanisms by which P can be induced 
by 00, therefore, we analyze within the BP approach 
a simple one-dimensional tight-binding model where 00 
can be controlled by a Hubbard-like interaction. We will 
consider only spinless fermions, which allow to easily dis- 
entangle the role of spin and charge from orbital degrees 
of freedom in causing a finite polarization. We will show 
that bond inequivalency without CO but with a finite 
orbital-occupancy disproportionation is enough to lead 
to a ferroelectric state, provided the intra-orbital hop- 
ping interactions are different. The choice of this ex- 
tremely simple model is motivated by the clear insight 
it provides within the BP framework, rather than by its 
straight application to some real material. Nonetheless, 
we expect the proposed mechanism to be rather general 
and possibly relevant for the recently proposed multifer- 
roic vanadium spinel CdV90zt[2^ and similar systems, as 
we discuss in our concluding remarks. 

Model and analytical formulas. The possible interplay 
between orbital degrees of freedom and bond inequiva- 
lency in inducing a bulk P can be addressed by consider- 
ing a spinless two-band Peierls-Hubbard model, defined 
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on a onc-dimensional chain, whose Hamiltonian reads 

K 



+ 4i^.i4^.2, (1) 



with dj^ creating an electron in the orbital 7 of site 
j. The hopping parameters are modulated by the 
e — I interaction, with coupling constant a, as tj''' = 
^77 [1 _[_ a{uj — Uj+i)] , where Uj are the displacement of 
the ions from their equilibrium position along the chain, 
with the associated clastic constant K. Without loss 
of generality, in the following we will neglect any inter- 
orbital hopping, assuming t'^'^ = 5^,^d^ with ti ^ t2. 
The last term describes the Coulomb interaction between 
electrons occupying different orbitals, parametrized by 
the Hubbard constant U . 

Let us consider the half-filling case with average lo- 
cal occupation fixed to n = 1[23|. Peierls-like dimer- 
ization can be realized via a staggered pattern of dis- 
tortions Uj = (— l)-'uo, with uo to be variationally de- 
termined, whereas the Coulomb interaction can be de- 
coupled via standard Hartree-Fock linearization. The 
model is then formally equivalent to a single-band spinful 
Pcicrls-Hubbard model, whose staggered magnetization 
m = {—^y{nj^ — nj^) maps to the difference between av- 
erage orbital occupancies A ~ (—1)'' {nji—nj2), that may 
serve as the order parameter for orbital ordering. Since 
the periodicity of dimerization and of staggered 00 is the 
same, the Hamiltonian in reciprocal space can be defined 
in the reduced Brillouin zone, reading 



(5^(fc) + (-l)^iC/A 



r(fc) 



(2) 



in each orbital sector, where e^{k) = — 2t^cosfc and 
(5^(fc) = i4Q; tip t-y sin fc, with both Uo,A to be self- 
consistently determined. Eigcnstates are given by: 



El{k) = \u± ^e,{kY + \5,{k)\^ + \ ([/A)^ (3) 

implying an insulating groundstate for nonvanishing mq 
or A, which at half filling must fulfill the following self- 
consistency equations: 
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^e,iky + \6,ikW + \{UA)^' 



(4) 
(5) 



Depending on the coupling constants a, U, four phases 
are found characterized by different structural and orbital 
properties. For small a and U — 0, we found a metallic 



ground-state, with no orbital-occupancy disproportiona- 
tion and no structural distortion. As U, a are increased, 
three different insulating ground-states develop with fi- 
nite A or wq. a first insulating phase, with A 7^ and 
Mo = 0, is induced by increasing the e — e interaction, ex- 
actly mapping onto the antiferromagnetic insulator found 
for the half-filled single-band Hubbard model. On the 
other hand, a dimerized state, with uq 7^ and A = 0, 
develops with increasing a at small U, which is the con- 
ventional Peierls insulator with structural distortions. A 
third insulating phase intrudes between the first two at 
moderate values of U by increasing the e — I interaction; 
this insulating state is characterized by the coexistence 
of 00 and structural dimerization. The phase diagram 
in the U ~ a space is shown in Fig. [1] for a given ra- 
tio between hopping parameters, namely ti/t2 = 3. The 
evolution of dimerization uq and orbital occupancy dis- 
proportionation A is also shown as a function of U, high- 
lighting the coexistence regions at different values of the 
e — I interaction. Even if the phase boundaries depend on 
the ratio ti/t2 and may be affected by the mean-field ap- 
proximation, the phase diagram is qualitative the same 
as the one obtained for the single-band Peierls-Hubbard 
model by more accurate numerical methods [2^. 

Turning to ferroclcctricity, the key quantity to be con- 
sidered is the overlap matrix defined as Sm,m'ik,k') = 
{'4'm,k\G~^^^\fpm' ,k'), which allows to evaluate the ex- 
pectation value of the position operator x on the eigcn- 
states ipm,k, with m the occupied band index, which obey 
periodic boundary conditions on a Irf chain of Icnght 
L|25j. Since the overlap matrix has nonvanishing ele- 
ments only when each pairs of k,k' differs by e = 2tt/L, 





FIG. 1: Left: Phase diagram in the U — a space for the two- 
band spinless Peierls-Hubbard with inequivalent hoppings 
ti/t2 = 3: the thick line represents the metallic solution 
found at [/ = 0, I and III label respectively the Peierls and 
00 insulator, with II the mixed-character insulator intrud- 
ing between them. All parameters are expressed in units 
of t = {ti + t2)/2 = 1 and the elastic energy is fixed to 
K — 1. Right: evolution of the dimerization parameter mo 
(thin lines) and of the orbital-occupancy disproportionation 
A (thick lines) as a function of U for different values of the 
e — I interaction. Regions where OO coexists with structural 
distortion are highlighted. 
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one can explicitly evaluate its elements on the eigcn- 
states of Hamiltonian ([2]) in the limit of L — oo, finding 
Sm,m' {k,k + e) ^ <5m,m' (1 + <S'm) with: 



i-ir 
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([/A)2 +4|J™(fc)| 



(6) 



where m = 1,2 labels the occupied bands with eigen- 
values E}^{k), E'^{k) as given in Eq. Q- By assuming 
an uniformly charged ionic background, no ionic contri- 
bution arises from the chain distortion in the absence of 
charge disproportionation; polarization, then, has only 
an electronic origin and is evaluated as[25j: 

-P = -Im log det S{k, k + e) ^ - ^ 5™ dk, 



(7) 



e being the electron charge. By inspecting Eqs. ([6]), ([7]) it 
comes clear that two conditions must be simultaneously 
fulfilled in order to have a nonvanishing P: i) both 00 
and structural dimerization, with uqA ^ 0, must occur, 
and ii) hopping integrals must be different in different 
orbital channels, i.e. ti ^ ^2, otherwise the two addends 
in the r.h.s of Eq. ([7]) would cancel out exactly leading to 
P = 0. The explicit expression for Sm{k), as given in Eq. 

also implies that the polarization can be reversed ei- 
ther by swapping the orbital occupancy between the two 
inequivalent sites in the doubled cell or by interchanging 
short and long bonds of the distorted chain. A paradig- 
matic evolution of P in the dimerized orbital-ordered 
state is shown in Fig. ([5]). Enforcing the correspon- 
dence between the two-band spinless and one-band spin- 
ful Peierls-Hubbard model, one immediately recognizes 
that the second condition is never fulfilled in the spin- 
ful model, since hopping integrals are spin independent 
and identical. This explains why the antifcrromagnetic 
ground state of the single-band Peierls-Hubbard model 
is never ferroelectric unless a CO is realized, e.g., by in- 
cluding a staggered potential, as already discussed in Rcf. 

A pictorial interpretation of the origin of the electronic 
P in the orbital-ordered dimerized state is provided by 
looking at the center of the Wannier functions (WF) of 
the two occupied bands, that are easily evaluated for the 
Id system under consideration. As pointed out in Ref. 
26l . in fact, the construction of maximally localized WFs 
is not needed in the present case; WF centers can be ob- 
tained as the eigenvalues of a matrix A, constructed as 
the product of the unitary parts of the S matrices along 
the k-point string (given by the matrix product VW^ 
taken from the singular value decomposition S = FEW^^, 
where V and W are unitary and S is a diagonal matrix 
with nonnegative diagonal elements) . In the limit of infi- 
nite chain (with e — >■ 0), it turns out that the WF centers 
rm are given by the two addends in the r.h.s of Eq. ([7]), a 
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FIG. 2: Evolution of the Berry-phase polarization (thick lines) 
as a function of U sd, a — 0.25 and ti/t2 = 3. Thin lines are 
the dimerization parameter uq (dot-dashed line) and the 00 
parameter A (solid or dashed); dashed lines correspond to the 
swapping of orbital occupancy by keeping fixed the structural 
distortion of the chain, causing a reversal of P. 



well-known result in the context of the modern theory of 
polarization [sj. That means that in the orbital-ordered 
dimerized state the two WF centers move in opposite 
directions along the chain, refiecting the staggered or- 
bital occupancy; the different hopping integrals ti ^ t2 
(implying \Si\ ^ \S2\), however, make the displacements 
of Tm with respect to their nonpolar positions inequiv- 
alent, giving rise to uncompensated local dipoles that 
sum up resulting in a bulk P. This is shown in Fig. |31 
where the evolution of is plotted as a function of U 
for a — 0.25 and ti/t2 ~ 3. In the Peicrls insulator, 
where bonds differentiate in alternating short and long 
ones, both the WF centers lie in the middle of the short 
bond. On the other hand, when 00 develops, e.g. with 
prevalent 7=1 character on odd sites and 7 = 2 on even 
sites, the corresponding WF centers move towards the 
most occupied sites with same orbital character. Due to 
the different values of hopping integrals, these displace- 
ments are inequivalent and cause the local dipoles to be 
uncompensated. Eventually, when the structural distor- 
tion disappears in the 00 insulator, each lies on a 
site of the chain, as schematically shown in Fig. O)). 

Conclusions. We analysed the ferroelectric properties 
of a simplified model in the framework of the modern the- 
ory of polarization in order to unveil the necessary condi- 
tions by which 00 phenomena may lead to an electronic 
bmP. The spinless two-band Peierls-Hubbard model al- 
lowed us to consider the cooperative interplay between 
00 and Peierls-like dimerization along a Id chain. The 
derived analytical formulas, even if obtained at a mean- 
field level, clarify that P may arise from 00 only if a con- 
comitant bond inequivalency (appearing as a structural 
dimerization in the considered model) develops. Naively, 
the emerging P may be interpreted as a superposition of 
two inequivalent charge-localization phenomena in dif- 
ferent orbital sectors. By looking at Wannier-function 
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plitudcs along a single one-dimcnsional "chain", e.g., 
parallel to yz direction would fulfill tyz.zx = and 
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FIG. 3: a) Evolution of the Wannier function centers rm, in 
units of lattice spacing ao = I, as a function of U for a = 0.25 
and t\/t2 = 3. b) Schematic picture of the displacement of 
WF centers as the system evolves from the Peierls insulator 
(I) to the OO insulator (III) through the dimerized orbital- 
ordered state (II), highlighting the inequivalent local dipoles 
that develops in phase II (empty circles indicate the undis- 
torted position of sites along the chain; the displacements are 
exaggerated to make the picture clearer). 



centers, that can be viewed as the localization centers of 
the continuous electronic charge distribution, local elec- 
tric dipoles develop as the WF centers, with given orbital 
character, move towards the site with most orbital char- 
acter of the same kind; in the presence of bond inequiv- 
alency, these dipoles are uncompensated, giving rise to 
a nonzero P. Finally, the polarization can be reversed 
by reversing 00 or the chain dimcrization. The unveiled 
mechanism is different from the one previously reported 
for undopcd manganites[l3l. In that case, Jahn- Teller 
lattice distortions tuned 00, whereas a large Hund cou- 
pling in the E-type antiferromagnetic phase constrained 
electrons to hop along ferromagnetic zig-zag chains; the 
change of electron motion around each site displaying 
00 was responsible, through direction-dependent inter- 
orhital hopping processes, for the phase change of the 
Bloch wavefunctions leading to P. 

Even though the analysis has so far been restricted 
to little more than a toy model, the proposed mech- 
anism may prove relevant for spinel crystals, such 
as MgTi204[23l, CuIr9 SA 12811 or vanadium spinels as 
ZnV9 0A[29{ or CdV2 04[23. In the simplest approxi- 
mation, the electronic structure of these systems has es- 
sentialy one-dimensional character, stemming from the 
strong direct d — d overlap of the i?— site t2g orbitals. 
The strong anisotropy of the d— electron distribution, in 
fact, is such that electrons in xy,yz,xz orbitals would 
hop preferably in spinels along xy,yz,xz directions re- 
spectively, implying also strongly direction-dependent 
(orbital-dependent) hopping integrals along each effec- 
tive one-dimensional "chain" [s^ . For instance, in vana- 
dates (V^"*") one finds that the xy— orbitals are always 
occupied due to a tetragonal distortion while the re- 
maining t2g electron moves in degenerate yz, zx bands; 
adopting Slatcr-Koster paramctrization, the hopping am- 



32| . close to our choice of pa- 
rameters. Each chain in the spinel structure may then 
be subject to Peierls transition, which has been sug- 
gested to be strongly tied to possible orbital ordering 
phenomenally 34 1, eventually driving the system insu- 
lator. Interestingly, at least one of the forementioned 
vanadium spinels, CdV2 04, has been reported to be 
ferroelectric [22[. Other issues should be taken into ac- 
count if willing to extend our analysis to more realistic 
situations, such as the forementioned spinels. In par- 
ticular, the role of spin degrees of freedom, that has 
been overlooked in order to highlight the role of orbital 
degrees of freedom, may play a relevant role in induc- 
ing the bond inequivalency, analogously to what hap- 
pens in CasCoMnOg or in collinear rare-earth mangan- 
ites i?Mn03. Furthermore, the interplay between spin 
and orbital degrees of freedom could be relevant for novel 
multiferroic or magnetoelectric effects. 
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